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Abstract. We analyse the sensitivity of IceCube-DeepCore to annihilation of neutralino dark matter 
in the solar core, generated within a 25 parameter version of the minimally supersymmetric standard 
model (MSSM-25). We explore the 25-dimensional parameter space using scanning methods based 
on importance sampling and using DarkSUSY 5.0.6 to calculate observables. Our scans produced 
a database of 6.02 million parameter space points with neutralino dark matter consistent with the 
relic density implied by WMAP 7-year data, as well as with accelerator searches. We performed a 
model exclusion analysis upon these points using the expected capabilities of the IceCube-DeepCore 
Neutrino Telescope. We show that IceCube-DeepCore will be sensitive to a number of models that 
are not accessible to direct detection experiments such as SIMPLE, COUPP and XENON-100, nor 
to current LHC searches. 
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1 Introduction 

Much of the development of supersymmetry (SUSY) was motivated by problems within the Standard 
Model (SM), such as the fine tuning of the Higgs mass [1, 2] and the unification of gauge coupling 
constants at high energy [3-7]. When the lightest neutralino is the lightest supersymmetric particle 
(LSP) and _R-parity is conserved, the lightest neutralino is a natural candidate for dark matter, 
providing a solution to one of the most important problems in astrophysics. 

Neutralinos can accumulate in the centre of astrophysical bodies such as the Sun, where they 
would annihilate into a range of SM particles. Subsequent decays and interactions of these SM particles 
would yield neutrinos, which have the unique ability to escape the solar core and reach detectors on 
Earth [8-10]. 

One such detector is the IceCube Neutrino Telescope [11]. Previous analyses of data recorded 
by IceCube in the 22- and 40-string configurations have already provided limits on neutralino dark 
matter annihilation in the Sun, the neutralino-proton cross-section, and the resulting muon flux[12, 
13]. Further papers have used projected capabilities of the full 86-string IceCube-DeepCore detector 
(IceCube-86) to analyse parameterisations of the MSSM in more detail, such as the four-parameter 
Constrained MSSM (CMSSM) [14, 15], the seven parameter MSSM-7 [16, 17], and the 19 parameter 
phenomenological MSSM (pMSSM) [18]. 

In this paper we present the results of an exploration of a 25 parameter version of the MSSM 
and the sensitivity of IceCube-86 to these models. This analysis involves two new elements compared 
to previous work: the parameter space has been enlarged to 25 parameters, reducing the artificial 
constrictions placed upon the MSSM; and the IceCube-86 model exclusion confidence level calculations 
now include point-by-point optimisation of the reconstructed cut angle around the Sun, which was 
not included in a similar analysis performed in [17]. 

The structure of this paper is as follows. In Section 2 we introduce the IceCube Neutrino 
Telescope and its salient characteristics. Then, in Section 3, we present the details of our analytical 
methods: the MSSM-25 parameterisation within which we worked (3.1), the method by which we 
explored this parameter space (3.2), the calculation of signal rates in IceCube (3.3), and finally 
the calculation of IceCube model exclusion confidence levels (3.4). In Section 4 the results of this 
exploration and analysis are presented and discussed. 

2 The IceCube Neutrino Telescope 

Completed in December 2010, the IceCube Neutrino Observatory consists of a cubic kilometre of 
extremely transparent natural ice deep below the South Pole, instrumented with 5160 detector optical 
modules (DOMs), along with 320 DOMs embedded in ice-filled tanks at the surface (known as the 



- 1 - 



IceTop array). The in-ice DOMs are arranged along 86 strings lowered down vertical wells melted 
into the ice using a hot water drill. Seventy-eight of the strings are arranged in a hexagonal grid with 
a horizontal spacing of 125 m between strings, with DOMs separated vertically by 17 m and situated 
between 1450 m to 2450 m below the surface. The remaining eight strings are clustered in the centre 
with a horizontal spacing of 72 m, and have 50 of their DOMs positioned with 7 m vertical spacing 
starting at the bottom. The remaining ten are arranged higher up the strings to act as a veto. These 
eight strings and the 12 adjacent standard IceCube strings make up the DeepCore subarray, and 
provide increased sensitivity at lower energies and reduce the energy threshold. The IceTop array 
is used to study cosmic ray flux and composition, and serves as a veto for high energy cosmic ray 
showers [11]. 

IceCube can reconstruct the type, energy and direction of the incoming neutrinos by observing the 
Cerenkov radiation emitted by leptons created in charged current interactions between the neutrinos 
and nuclei in the ice. Muons produced in such interactions are of particular interest, as their mean 
free paths can be up to 10km for the most energetic neutrinos. One can thus detect neutrinos that 
interact outside the volume of instrumented ice making up IceCube, and so the effective volume for 
muon neutrinos exceeds 1km 3 [11]. The IceTop array is used to study cosmic ray flux and composition, 
and serves as a veto for high energy cosmic ray showers [11]. 

The sensitivity of IceCube to neutrinos of different energies is characterised by the effective 
area. For this study we use the effective area for muon neutrino detection by the 86-string IceCube 
configuration (IceCube-86), derived from simulation and presented in [19]. The sensitivity study 
presented in [19] is for the full 86 string detector with the initially proposed DeepCore geometry, that 
consisted of only six additional strings instead of the eight deployed. The study was performed as 
a full analysis in all details on simulated backgrounds, including realistic data processing and event 
selection. After the work presented in this paper was completed, preliminary results from the IceCube- 
79 detector were released [20]. A comparison of the effective area from [19] and these preliminary 
results suggests a moderate over-estimation of the effective area by the earlier sensitivity study: a 
factor of ~3 at 100 GeV neutrino energy, dropping to ~1.5 at 1 TeV. Given the order unity difference, 
the effective area from the IceCube-86 sensitivity analysis [19] will suffice for our purposes here. 

3 Methods 

3.1 The MSSM-25 Parameterization 

The full MSSM has 178 free parameters: 158 from the soft SUSY-breaking sector, one from the cou- 
pling of the two Higgs supermultiplets, and the remaining 19 from the Standard Model [21]. 54 of 
these parameters can be removed by suppressing the C-terms of the soft SUSY-breaking Lagrangian, 
leaving 124 parameters. We direct the reader to [21] for more information. For phenomenological 
studies it is useful to reduce this number in order to make exploration of the parameter space compu- 
tationally feasible. This reduction is performed via a series of assumptions that become progressively 
more severe as more parameters are eliminated. Previous analyses of neutralino dark matter signals 
in IceCube have used models with four, seven, nine and 19 parameters [14-18]. For this research we 
used a 25-parameter version of the MSSM. 

The potential utility in increasing the number of parameters lies in reducing the artificial constric- 
tions placed upon the MSSM and allowing different combinations of parameters that could perhaps 
yield dark matter candidates with previously unseen properties. The disadvantage of increasing the 
number of parameters is the increased computation time needed to sample the parameter space. 

The parameters of the MSSM-25 are listed in the left column of Table 1. The Higgs sector is 
described by the ratio of Higgs VEVs tan j3, the CP-odd Higgs mass itia generated after electroweak 
symmetry breaking, and the /i parameter from the MSSM superpotential. The gaugino sector of the 
soft SUSY-breaking Lagrangian £ so ft is parameterised by the three mass terms Mi, M2, and M3. 

The five Hcrmitian 3x3 sfermion mass-squared matrices in £ so ft are simplified by assuming the 
absence of flavour changing neutral currents (FCNCs), forcing all non-diagonal entries to zero. Each 
mass squared matrix therefore only has three real free parameters, and so the sfermion sector has a 
total of 15 parameters. 
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For the triple scalar couplings in £ so ft we again apply the prohibition on FCNCs to eliminate 
the off-diagonal entries. We also restrict the matrices to be Hermitian (and so the diagonal elements 
to be real) in order to eliminate any CP-violating phases [22]. We also make the assumption that the 
elements of the triple scalar coupling matrices are proportional to the corresponding elements of the 
Yukawa matrices y u , yd, and y e . We can set the first and second entries of a u and ad to zero as the 
corresponding Yukawa couplings are negligibly small, but we retain one parameter, a ol / 2 , for the first 
and second entries of a e , as it is relevant to the calculation of the anomalous magnetic moment of the 
muon [23]. We finally arrive at the following triple scalar coupling matrices: 



\ / 

a„ = i a d = 

OnaKas / V 



3.2 Parameter Scans 





We performed scans of this parameter space using importance sampling with ADSCAN [24], an 
adaptive scanning programme based on the VEGAS algorithm [25]. The VEGAS algorithm works 
by performing an initial random scan of the parameter space and then focusing subsequent scans 
on areas that produce important values, as defined by an importance function G. Our importance 
function was 



G (fl x h 2 ) = exp 



1 f fi^/l 2 — 



(3.2) 



where = 0.1120, the cold dark matter relic density derived from 7-year WMAP measure- 

ments, with lcr uncertainty of ±0.0056 [26]. For the error in the denominator we chose a^h 2 = 0.01 
to allow for some theoretical error above the lcr experimental error of ±0.0056, as we did in the 
corresponding MSSM-7 ADSCAN parameter scans in [17]. The parameter ctq^ determines the width 
of the Gaussian importance function used in the scanning process, with larger numbers generating a 
larger scan. 

We ran the adaptive scanning programme multiple times with varying parameter limits, which 
are listed in Table 1. We rejected points in the parameter space that gave negative mass-squared 
terms, and where the lightest neutralino was not the Lightest Supersymmetric Particle (LSP). We 
eliminated additional points by applying experimental limits on the accelerator observables produced 
by each point: the sparticle mass spectrum [27]; standard model and SUSY Higgs masses from LEP 
and CDF measurements [27]; the b — > 57 branching ratio using results from the Heavy Flavour 
Averaging Group [28]; the invisible Z° boson decay width [27]; and the magnetic moment of the 
muon derived from data from the BaBar and KLOE experiments [29-31]. 

As the primary focus of this paper is the impact of IceCube-DeepCore, our overall strategy was to 
deliberately lean towards a more inclusive scan, choosing to discard only points that disagreed with 
the most established and well-understood experimental limits. We made this choice so as to keep 
scans simple and easy to understand, as this helps to more clearly focus on the specific impacts of 
IceCube-Deepcore on the MSSM-25. In this way we avoid worrying about erroneously discarding some 
models due to e.g. the unknown model-dependence of existing LHC SUSY limits when translated into 
the MSSM-25, or the massive amount of computational resources required to obtain this knowledge. 

From the database of points satisfying physicality, the demand for a neutralino LSP, and accel- 
erator bounds, we selected those that gave a relic density within the 2cr error bounds of the 7-year 
WMAP measurement, i.e. in the range 0.1120 - 0.0112 < Vt x h 2 < 0.1120 + 0.0112. We found a total 
of 6.02 million such points. 

For calculating the values of all observables, we used the DarkSUSY 5.0.6 software package 
[17, 32], with default options except where otherwise specified. In particular, for the calculation of 
neutralino-nucleon scattering cross-sections this implies A s = —0.12 for the strange quark content of 
the nucleus, and S^at = 45 MeV, cr = 35 MeV for the hadronic matrix elements. 
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Parameter Units Limits 



Gaugino mass terms M 1; M 2 , and M 3 


GeV 


-16000 


16000 


Higgs parameter /i 
CP-odd Higgs mass tua 
Higgs VEV ratio tan j3 


GeV 
GeV 


-100000 
-8000 
2 


100000 
8000 
65 


Sfermion mass parameters , , 

™q 3 > m iv m= 2 , mg 3 , m^, m^, m & , 

m L^ m L 2 ' m L 3 ' TO 5i> TO 5 2 ' and m s 3 


GeV 


100 


25000 


Triple scalar couplings <z U 3, ad3, <2 e i/2> and a e 3 


GeV 


-50000 


50000 



Table 1. MSSM-25 Parameter Space Maximum Exploration Limits 



3.3 IceCube Signal Rates 

We post-processed the 6.02 million points with suitable relic density using DarkSUSY 5.0.6, as de- 
scribed in [17], to find the IceCube signal rates that would be produced by each model. 

We start by modelling the accumulation and annihilation of neutralinos in the Sun. The total 
population of neutralinos in the Sun N(t) is described by the equation [33] 

^ =C c -C a N{t)\ (3.3) 
di 

where C c is the constant neutralino capture rate and C a is defined as 

C a = {a a v)^ with V 3 ■« 6.5 x 10 28 C^^) ^cm 3 . (3.4) 

The Vj terms are effective volumes, which take into account the quasi-thermal distribution of neu- 
tralinos in the Sun. 

The capture rate calculation depends on the halo density and velocity profile of dark matter, the 
neutralino mass and interaction cross-section, and the properties of the Sun; see e.g. [34] for a detailed 
description. We assume a standard dark matter halo model, with the Sun moving at vq = 220kms _1 
through a halo with local dark matter density po — 0.3 GeVcm~ 3 , and dark matter velocities following 
a Maxwell-Boltzmann distribution with dispersion v — 270 km s . We did not include the detailed 
effects of diffusion and planets upon the capture rate, as the simple free-space approximation [35] 
included in DarkSUSY turns out to be highly accurate [36]. The annihilation rate of neutralino pairs 
is given by 

r o (t) = \c a N 2 (t). (3.5) 
Solving Equation 3.3 gives us the annihilation rate at a given time t as 

r o (i) = ^tanh 2 ('^ (3.6) 

where 

1 



\J C c C a 



(3.7) 



is the capture-annihilation equilibrium time scale. Taking t as the age of the solar system (4.5 billion 
years) and approximating the structure of the Sun to be constant over its lifetime gives the current 
neutralino annihilation rate in the Sun. 



- 4 - 



From the annihilation rate T a (t) and branching fractions, DarkSUSY calculates a predicted neu- 
trino flux spectrum at the IceCube detector, using lookup tables previously computed with WimpSim 
[37]. WimpSim takes neutrino production in the Sun and propagation through the Sun, Earth and 
interplanetary space into account using a full three-flavour Monte Carlo simulation. Together with 
the effective area of IceCube-86 and the mean angular error obtained from [19], and some choice of 
analysis region around the Sun, the calculated flux yields an expected signal rate for any given model. 

The other necessary ingredient for our model exclusion calculations is the number of background 
events. We generated these using a bootstrap Monte-Carlo re-simulation of the expected background 
rates away from the Sun, as derived in the sensitivity study [19]. The total expected number of 
background events over the whole sky was 15552 [38], consisting of atmospheric neutrino and muon 
backgrounds. Given only the total number of background events, we assumed all events to be uni- 
formly distributed in declination. This is a good approximation for atmospheric neutrinos within the 
narrow declination band of the analysis, and a conservative choice for the down-going atmospheric 
muon component. This results in an overall robust background expectation. The single realisation 
of this expected background that we employed for this study contained 1452 simulated events within 
20° of the solar centre. 

3.4 Likelihood Calculation 

In this work we perform a model exclusion analysis similar to that of Section 4.3 in [17]. Compared 
to [17], we employ an expanded parameter space and optimise our angular event cut around the solar 
position. We also look in far more detail than [17] into the impacts of IceCube-86 on the model 
parameters and their corresponding derived observables. To quantify the exclusion CL we calculated 
the p- value for each model using the method outlined in Section 3.6 of [17], specifically Equation 3.28 
and its contributors. The CL at which a model can be excluded is then 1 — p. 

We excluded all events with reconstructed angle outside some bin of width <j> centred on the 
Sun. We optimised the cut angle by repeating the p- value calculation for <j) = 3°, 6°, 9°, and 20°, and 
selecting the lowest p-value (i.e. highest exclusion CL) for each model. We optimised the cut angle 
because different models produce signals with different angular distributions. For instance, higher 
mass neutralinos tend to produce higher-energy neutrinos, which are more accurately reconstructed 
by IceCube and are therefore more densely clustered around the solar position. Using a small cut 
angle in this case reduces the number of background events included in the bin more than it reduces 
the number of signal events. A brief analysis of optimal cut angle is presented in Section 4. 

Because we work entirely with simulated data in this paper, we can simply optimise the cut 
angle by choosing the value that produces the smallest p-value. Working with real data however, 
selecting the cut angle that gave the lowest p-value would constitute an a posteriori analysis choice, 
which would make the statistical treatment less than rigorous. As we show later, the cut angle giving 
the lowest p-value is strongly correlated with the neutralino mass and somewhat correlated with the 
annihilation channel, the quantities that dchne the expected neutrino energy spectrum (this is to be 
expected, because the angular distribution of muons is heavily dependent upon the neutrino energies). 
In an analysis of real data, one should therefore optimise the cut angle in advance using the expected 
neutrino energy spectrum for a given model, and an ensemble of background-only simulations like the 
one we employ here. 

4 Results and Discussion 

Figure 1 shows the results of the model exclusion analysis we performed on our set of 6.02 million 
points. Care must be taken in interpreting these plots, as points of higher exclusion CL are plotted 
on top of points with lower exclusion CL, and so the former can obscure the latter. An area of cyan 
or blue on Figure 1 means that IceCube has exclusion (or detection) capability of 3a (99.7% CL) 
or better, for a certain range of interaction cross-sections, neutralino masses, and other MSSM-25 
parameters. Nevertheless, one can see in the areas of blue and cyan points that IceCube has the 
ability to exclude models at better than 99% CL well beyond even the 90% CL reach of current direct 
detection experiments such as SIMPLE [39], COUPP [40] and XENON-100 [41, 42]. 
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Lightest neutralino mass ro, (GeV) Lightest neutralino mass m x (GeV) 



Figure 1. Spin-dependent neutralino- proton cross-section asu, P (left), and spin- independent neutralino- 
proton cross-section <tsi,ji (right) as functions of lightest neutralino mass m x , for points derived from explo- 
rations of the MSSM-25 parameter space. In the left panel 90% CL Spin-dependent WIMP-proton cross-section 
limits from SIMPLE [39] and COUPP [40] direct detection experiment are displayed as magenta and black 
lines respectively. In the right panel the 90% CL Spin-independent WIMP-nucleon cross-section limit from 
225 live days of XENON100 direct detection experiment data is displayed as black line [41]. Colour coding 
indicates predicted IceCube-86 model exclusion CL. The areas of cyan and blue points show that IceCube-86 
has the ability to exclude models beyond the reach of current direct detection experiments such as SIMPLE, 
COUPP, and XENON100. 




annihilations (km 2 year l ) annihilations (km 2 year ') 

Figure 2. Spin-dependent neutralino-proton cross-section asu, P (left) and spin-independent neutralino-proton 
cross-section asi, P (right) against muon flux in IceCube from neutralino annihilations, for points derived from 
explorations of the MSSM-25 parameter space. Colour coding indicates predicted IceCube-86 model exclusion 
CL. 

In Figure 2 we see spin-dependent and spin-independent neutralino-proton cross-sections plotted 
against the predicted muon flux in IceCube-86 from solar neutralino annihilations. For these calcula- 
tions we used a muon energy threshold of 1 GeV, and a maximum angular separation of 30° between 
the solar centre and moun arrival angle. 

In general models which can be excluded at higher CL have higher muon fluxes - high muon flux 
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Muon flux in IceCube from neutralino 
annihilations (km -2 year -1 } 

Figure 3. Lightest neutralino mass m x against muon 
flux in IceCube from solar neutralino annihilations, 
for points derived from explorations of the MSSM-25 
parameter space. Colour coding indicates predicted 
IceCube-86 model exclusion CL. The effect of the en- 
ergy dependence of IceCube's effective area can be 
seen in the distribution of high exclusion CL points: 
lower mass neutralinos will yield lower energy neutri- 
nos/muons, and so a higher flux is necessary to create 
the same signal and thus exclusion CL in IceCube. 




SD neutralino-proton cross-section a SDp (cm 2 ) 

Figure 4. Spin-independent neutralino-proton cross- 
section o"si,p against spin-dependent neutralino-proton 
cross-section aso, P for points derived from explorations 
of the MSSM-25 parameter space. Colour coding 
indicates predicted IceCube-86 model exclusion CL. 
Dashed black line indicates equality between spin- 
independent and spin-dependent neutralino-proton 
cross-sections. 



leads to high signal, and high signal models are easier to rule out. This relationship, however, is not 
exact, as witnessed by the fact that the different regions of exclusion CL (i.e. different coloured points) 
are not cleanly banded. In Figure 2 one can see points with lower exclusion CL that nonetheless have 
higher muon fluxes than points with higher exclusion CL. This is a result of IceCube's effective area 
increasing with energy, as seen in Figure 2 of [19]. A model can have a high muon flux but if it also 
has a comparatively low neutralino mass, then it will produce lower energy neutrinos and muons, 
which leads to a lower IceCube signal rate. Figure 3 throws this interplay into stark relief; models 
with lightest neutralinos of lower mass will produce neutrinos with lower energy, and so compared to 
models with higher lightest neutralino masses a higher flux is necessary to compensate for the reduced 
effective area at lower energies. 

Figure 4 compares spin-independent to spin-dependent neutralino-proton cross-sections. Note 
that many models have spin-dependent neutralino-proton cross-sections that are larger than their 
spin-independent neutralino-proton cross-sections. The majority of models with high exclusion CL 
have high spin-dependent cross-sections. The long vertical band of high exclusion CL at high spin- 
dependent cross-section are regions where spin-dependent scattering is the dominant capture mode of 
neutralinos in the sun. As one moves upwards along this band a 'turn-off' is reached at approximately 
a si.p — 10~ 44 cm 2 ; above this point, spin-independent nuclear scattering is a significant source of 
neutralino capture in the Sun, and even dominates the capture rates in many cases. 

Figures 5 to 8 show the exclusion CL across a selection of parameters. Several parameters display 
a strong influence on exclusion CL, as witnessed by concentrations of blue and cyan points. These 
concentrations can be found where 100 GeV < Mi < 200 GeV, (Figure 5), around /i = 200 GeV 
(Figure 7 Left and Right), and in the region bounded roughly by 100 GeV < < 500 GeV and 
400 GeV < n < 1200 GeV (Figure 7 Left). There is a clustering of high exclusion CL points at high 
values of the triple scalar couplings, as seen in Figure 6. The sfermion mass parameters and tan /3 
show no influence over exclusion CL, as can be seen from the roughly even distribution of blue and 
cyan points over their entire ranges. 

Here it is useful to recall our earlier note at the start of Section 4 regarding the interpretation of 
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Figure 5. Wino mass parameter M2 against bino mass 
parameter M\ for points derived from explorations of 
the MSSM-25 parameter space. Colour coding indi- 
cates predicted IceCube-86 model exclusion CL. Simi- 
lar results are given by M\ vs M3. 



Figure 6. Triple Scalar Couplings Od3 against a„3 
for points derived from explorations of the MSSM-25 
parameter space. Colour coding indicates predicted 
IceCube-86 model exclusion CL. Similar results are 
given by other triple scalar couplings. 
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Higgs Parameter ^ (GeV) Higgs Parameter ^ (GeV) 



Figure 7. CP-odd Higgs mass wia (left) and ratio of Higgs VEVs tan j3 (right) against Higgs parameter 
(i for points derived from explorations of the MSSM-25 parameter space. Colour coding indicates predicted 
IceCube-86 model exclusion CL. 



these plots. Points of high exclusion CL are plotted over those of lower exclusion CL, i.e. below the 
clusters of blue and cyan points there will be red and green points. Some parameter values produce 
large numbers of models that can be excluded at high CL by IceCube-86, but one cannot claim that 
IceCube-86 can exclude all points with those parameter values at high CL. The only correct way to 
make a statement about which parts of a parameter space are preferred or excluded, and at what CL, 
is to perform a full statistical global fit, as described in [17]. 

The distribution of points across the parameter space shows some distinct edges at certain 
parameter values: at M2 = 4000 GeV in Figure 5, at uia = 4000 GeV in the left panel of Figure 7, 
and the 'box within a box' features seen in Figures 6 and 8. These are artefacts of the scanning method 
and the parameter limits for various scans. For instance, we performed scans with limits on M 2 of 
±4000 GeV and ±16000 GeV, and so there are sharp cut-offs at these values. Similarly, we performed 
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10 J 10" 
1st gen. L-type squark mass parameter (GeV) 



10' 10 J 10" 

1st gen. L-type slepton mass parameter mi, (GeV) 



Figure 8. Second generation L-type squark mass parameters Mq 2 against first generation L-type squark mass 
parameter Mq (left), and second generation L-type slepton mass parameter against first generation L- 
type slepton mass parameter (right) for points derived from explorations of the MSSM-25 parameter 

space. Colour coding indicates predicted IceCube-86 model exclusion CL. Similar results are given by other 
squark and slepton mass parameters respectively. 




Lightest neutralino mass m x (GeV) SD neutralino-proton cross-section (j SDp (cm 2 ) 

Figure 9. Gaugino fraction ratio Z g /(1 — Z g ) of the lightest neutralino against lightest neutralino mass m x 
(left) and spin-dependent neutralino-proton cross-section asD.p (right), for points derived from explorations 
of the MSSM-25 parameter space. Colour coding indicates predicted IceCube-86 model exclusion CL. Note 
that the largest spin-dependent cross-sections, and therefore the best IceCube exclusion CL, occur where the 
lightest neutralino is a mixed gaugino-Higgsino. 

scans with limits on rtiA of ±4000 GeV and ±8000 GeV, leading to similar features in Figure 7 Left. 
We performed a large number of scans with sfermion mass parameters limited to ±16000 GeV, but 
only a few with limits expanded to ±25000 GeV. Thus the region bounded by this first limit was 
much more heavily sampled, and so has a higher density of points. Similarly the majority of scans we 
performed had limits on triple scalar couplings a U 3, ad3, a c i/2i an d a C 3 set to ±32000 GeV, but only 
a few had these limits expanded to ±50000 GeV. 

As shown in the left panel of Figure 9, IceCube has its strongest exclusion capability in the region 
where the neutralino is approximately equal parts gaugino and Higgsino. This corresponds to high 
spin-dependent neutralino-proton cross-section, as seen in the right panel of Figure 9. This agrees 




Figure 10. Lightest squark mass (left) and lightest neutralino mass (right) against giuino mass for points 
derived from explorations of the MSSM-25 parameter space. In the left panel the points to the bottom left 
of the magenta line are excluded by ATLAS at 95% CL, based on searches for direct production of coloured 
sparticles and their decay into jets and missing transverse energy (MET). The limits are derived by assuming 
that all sparticles are heavy (m = 5 TeV) except for the giuino and first- and second-generation squarks (the 
masses of which are scanned over), and the lightest neutralino (which is taken as approximately massless) 
[43] . Note that our MSSM-25 models do not rely on this assumption. In the right panel, points to the bottom 
left of the magenta line are excluded by CMS at 95% CL. This line is the combination of the best simplified 
model spectrum limits from the 'razor' [44] and Mt2 [45] analyses on giuino pair production and decay to 
diquarks + MET. These simplified analyses assume that each giuino in the pair decays identically directly to 
qqx°, i.e. that there are no intermediate states causing cascade decays of gluinos to neutralinos. This implies 
that squarks, charginos and the other neutralinos are either heavier than the giuino, or almost degenerate 
in mass with the lightest neutralino. Note that our MSSM-25 models do not rely on any such assumptions. 
Colour coding indicates predicted IceCube-86 model exclusion CL. 
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Figure 11. Optimal <f> cut against spectral hardness (left) and lightest neutralino mass (right) for points 
derived from explorations of the MSSM-25 parameter space. Colour coding indicates predicted IceCube-86 
model exclusion CL. Lighter and softer models typically benefit from larger cut angles. 



with the findings found in Section 4.3 of [17], where a similar analysis was performed using MSSM-7 
models. 

Figure 10 Left shows the lightest squark mass against the giuino mass for the points found during 
the MSSM-25 parameter space exploration. Also shown is the 95% CL exclusion limit derived from 
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4.71 fb _1 of data at = 7TcV from the ATLAS experiment at the Large Hadron Collider (LHC) 
[43]. This limit assumes a simplified set of SUSY models where all sparticles are given masses of 
5 TeV except for the gluino, first- and second-generation squarks, and the lightest neutralino. The 
first- and second-generation squark masses are scanned over, and the lightest neutralino is taken to be 
approximately massless. The ATLAS analysis also assumes that the coloured sparticles are produced 
directly and have direct decays into jets and missing transverse energy. Points below and to the left 
of the ATLAS limit curve are excluded under these assumptions. The MSSM-25 models we have 
produced are not based upon any such assumptions. The large number of blue and cyan points above 
and to the right of this limit curve show that IceCube-86 is sensitive at better than 99% CL to many 
models beyond the 95% CL exclusion capability of ATLAS. 

Figure 10 Right shows the lightest neutralino mass against the gluino mass for our database of 
MSSM-25 points. Also shown is the 95% CL exclusion limit derived from 4.7 fb -1 of data taken at 
■s/s = 7 TeV by the CMS experiment at the LHC, with points below this line considered excluded 
[44, 45]. The line is the combination of the simplified model spectrum limits from the 'razor' [44] 
and M T2 [45] analyses on gluino pair production and decay to diquarks + missing transverse energy 
(MET). The former is based on limits from 4 6-tagged jets -I- jets + MET shown in Figure 6 Top Left 
of [44], and the latter from the g — > bbx° channel using the M^b analysis, shown in Figure 6 Top 
Right of [45] . Such simplified model spectra assume that the gluinos in the pair decay identically, and 
that their decays to neutralinos proceed directly via g — ¥ qqx° only, without any cascade decays via 
intermediate states. This implies that other sparticles are either heavier than gluinos, uncoloured, 
or almost degenerate in mass with the LSP (such that their decay to neutralinos is kinematically 
suppressed or forbidden) . In practice this makes such limits applicable to the MSSM-25 only in cases 
where each squark, neutralino and chargino is either heavier than the gluino, or very close to the LSP 
mass; our MSSM-25 models do not make this assumption, so the CMS limit must be interpreted with 
quite some care. Nonetheless, the large number of blue and cyan points above and to the right of this 
limit curve suggest that IceCube-86 has the capability to exclude many models at better than 99% 
CL that are beyond the current 95% CL reach of CMS. 

Of the 6.02 million points we analysed, 5.75 million had an optimal cu t =3°, 32 had </> cu t = 6°, 
156,493 had cut = 9°, and 117,713 had <^ cut = 20°. Figure 11 Left shows the distribution of optimal 
(f> cu t for different values of the 'spectral hardness', defined as 

BF(t+t- , ti, W+W~ , ZZ) + ±BF(ZH,W ± H^) , 

Spectral Hardness = ; = = ■ ■ - — : , — — ; -. (4.1) 

BF(uu, dd, cc, ss, bb, e+e~, fi+fi~, HH) + \BF(ZH, W±H^) K ' 

Models with optimal 4> cu t = 3° spanned the widest range of hardness, but models with optimal 
<fi cut — 20° span a smaller range biased towards softer, lower energy spectra, as we might expect given 
that the angular resolution of IceCube is better at higher energies. 

Figure 11 Right shows the distribution of optimal <p cu t for different neutralino masses. Models 
with optimal cu t = 3° spanned the widest range of neutralino masses, while models with cut = 
20° spanned only a smaller range at lower neutralino mass. This confirms our earlier statement in 
Subsection 3.4: higher-mass neutralinos tend to produce higher energy neutrinos, which are more 
accurately reconstructed by IceCube and therefore more densely clustered around the solar position. 
In this case, a smaller cut angle reduces the number of background events more than it reduces 
the number of signal events. Conversely, lower-mass neutralinos will tend to produce lower energy 
neutrinos and so have a wider angular dispersion of reconstructed events; a wider cut angle is then 
required to maximise the number of signal events included in the analysis. 



5 Conclusions 

Our explorations of the MSSM-25 parameter space have yielded 6.02 million parameter space points 
that produce neutralino LSPs, satisfy constraints imposed by accelerator data, and produce a neu- 
tralino relic density within the 2a bounds of the 7-year WMAP measurement. Calculations of model 
exclusion CL for these points show that the 86-string configuration of IceCube will be able to rule 
out or detect at 3er (99.7% CL) or better a significant number of models that are beyond the 90% CL 



- 11 - 



exclusion limits of current direct detection experiments such as SIMPLE, COUPP, and XENON-100, 
and beyond the current 95% CL exclusion limits of ATLAS and CMS. 
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